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Abstract 

We use a recently proposed diagonalization/Monte Carlo computa- 
tional scheme to study relativistic two-body and three-body bound states 
in ((f) 6 — 4 )i+i theory. We find that the approach is well-suited for 
calculating bound state energies and wavefunctions. 

1 Introduction 

The structure and phenomenology of relativistic bound states are of primary 
interest to particle and nuclear physics. Unfortunately even in the simplest 
field theory models our understanding of relativistic bound states is limited 
by a number of complications. Monte Carlo simulation on the lattice is the 
standard tool for determining particle masses. One can find the low energy 
spectrum from the decay constants of Euclidean correlation functions. Though 
the procedure is generally successful it becomes more difficult and less reliable 
if the state of interest has the same quantum numbers as another lower energy 
state. Complications also arise for non-confining interactions, where states just 
above the continuum threshold can obscure the signal of the bound state. Both 
of these problems arise even in seemingly simple systems such as we consider 
here, namely, two-body and three-body bound states in ((f> 6 — </> 4 )i + i. 

One remaining complication is of a more operational nature, and it concerns 
the interface between computational methods and phenomenology. Much of 
the interesting and important physics associated with bound states is contained 
in subtle details. A measurement of gross properties such as the binding energy 
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is a necessary first step, but one also hopes to gain more descriptive information 
like the structure of the wavefunction, its response to various interactions, or 
point-by-point comparisons with intuitive non-relativistic models. 

These issues were part of the motivation for a recently proposed approach to 
computational quantum problems which combines diagonalization and Monte 
Carlo techniques jl], ^ . Diagonalization makes it possible to consider systems 
with sign oscillations and extract detailed information about wavefunctions and 
excited states, while Monte Carlo allows one to handle the exponential increase 
in the number of basis states for large volume systems. The first half of the 
method involves finding and diagonalizing the Hamiltonian restricted to an op- 
timal subspace. This subspace is designed to include the most important basis 
vectors of the lowest energy eigenstates. Once the most important basis vectors 
are found and their interactions treated exactly, Monte Carlo is used to sample 
the contribution of the remaining basis vectors. In this paper we demonstrate 
how the new diagonalization/Monte Carlo scheme is applied to the study of 
relativistic bound states. 

The essential problem of the lattice Monte Carlo approach to bound states 
near the continuum threshold involves finding the inverse Laplace transform of 
two-point correlation functions at energies close to thresholds where the spectral 
density is large. This is usually done by a non-linear multi-parameter fit of two- 
point correlators for several different operators. While in principle one might 
reach arbitrary accuracy by adding more parameters to the fit, in practice the 
presence of even very small stochastic error will destabilize the fit and imposes 
an upper limit on the accuracy. The main advantage of our method, therefore, 
is not in computing speed but in numerical stability. Our method also provides 
additional information about the wavefunction of the bound state. These two 
statements are not unrelated. Our method has greater numerical stability be- 
cause we are implicitly using the information contained in the wavefunctions 
to separate states. It is also worthwhile noting that in two dimensions the di- 
agonalization/Monte Carlo approach can be applied within the framework of 
discretized light cone and might prove useful since, aside from zero-mode prob- 
lems, the exponential growth of the number of basis states for large volume 
systems has so far provided a limit to discretized light cone. 

The organization of the paper is as follows. We begin with a discussion of 
the </> 6 — 4> Hamiltonian in 1 + 1 dimensions and its decomposition into mo- 
mentum modes. We then review the basic features of quasi-sparse eigenvector 
diagonalization and stochastic error correction. In the following two sections 
we briefly set aside the numerical discussion and derive solutions to the non- 
relativistic versions of the two-body and three-body bound state problems. We 
then compare our numerical results for the relativistic bound states energies 
and wavefunctions with an assortment of different approximations and results 
from the literature. 



2 



2 Fock states and the Hamiltonian 



Our starting point is the Hamiltonian density for <f> — <f> theory in 1 + 1 dimen- 
sions, 

where normal ordering is with respect to the mass fi. This theory was intro- 
duced in j| as a simple model with a stable bound state. In 1 + 1 dimensions, 
normal ordering is sufficient to renormalize the theory. 

We will consider the theory in a periodic box of length 2L. This allows 
us to expand in free particle momentum modes and rewrite the Hamiltonian in 
terms of mode creation and annihilation operators [[l], ^) , 

n 7i\-\ (-ri4— 



A ■ ( a "i+ a -n 1 ) K 6 +aL„ 6 ) 

32L 2 



raiH h«6=0 



In (Eh we have used 



Wn(A*) = (3) 

We can represent any momentum-space Fock state as a string of occupation 
numbers, \o-pf max , • • • , ojv max ), where 

a ri a n |0-AT max , ' ' ' , OAT max ) = O n |0-JV majt , ' ' ' , OAT max ) . (4) 

For the calculations presented here, we set the length of the box to size L — 
2.57r/i _1 and restrict our attention to momentum modes such that |n| < N mSLX , 
where iV max = 10. This corresponds with a momentum cutoff scale of A = 4/i. 
A change in the cutoff A^ max will lead to adjustments in /i and A in order to 
keep the physics independent of A^ max . However, these corrections begin at 
C(l/A^ lax ) as can be checked by power counting. The dependence on A^ max has 
been measured explicitly for the bound state data presented in this paper. It 
is considerably smaller than the dominant error, which is due to higher order 
corrections in the series stochastic error method, see below. 



3 QSE diagonalization 

In order to find approximate eigenvalues and eigenvectors of the Hamiltonian 
in (H), we use a method called quasi-sparse eigenvector (QSE) diagonalization. 
QSE diagonalization is an iterative algorithm that finds the most important 
basis vectors of the low energy eigenstates for a quantum Hamiltonian. We 
follow closely the discussion presented in 0] . We start with a complete basis for 
which the Hamiltonian matrix is sparse. In our case we use the momentum- 
space Fock basis. The QSE method involves the following iterated steps: 
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1. Select a subset of basis vectors {e^, • • • ,ei n } and call the corresponding 
subspace S. 

2. Diagonalize H restricted to S and find one eigenvector v. 

3. Sort the basis components of v according to their magnitude and remove 
the least important basis vectors. 

4. Replace the discarded basis vectors by new basis vectors. These are 
selected at random from a pool of candidate basis vectors which are con- 
nected to the old basis vectors through non- vanishing matrix elements of 
H. 

5. Redefine S as the subspace spanned by the updated set of basis vectors 
and repeat steps 2 through 5. 

If the subset of basis vectors is sufficiently large, the exact low energy eigen- 
vectors will be stable fixed points of the update process. This is the single 
eigenvector version of the algorithm. The multiple eigenvector version simply 
involves selecting more than one eigenvector each iteration and retaining the 
most important basis vectors for each of the eigenvectors. One of the most at- 
tractive features of our method is that we are working directly with the infinite 
dimensional Hilbert space. There is no restriction on the maximum occupation 
for any mode. This is a clear advantage over standard diagonalization methods 
such as Lanczos. 

In this analysis we work in the rest frame of the bound state and consider 
Fock states with total momentum zero. Due to <f> — ► —<f> reflection symmetry 
we can diagonalize the even and odd 4> sectors separately. We use the multiple 
eigenvector algorithm to obtain the four lowest eigenstates in the even sector 
and in the odd sector]^] For each QSE iteration, we select 50 new basis vectors 
and keep 400 basis vectors. The final QSE results are used as the starting point 
for the stochastic error correction method. 



4 Stochastic error correction 

The QSE method approximates the low energy eigenvalues and eigenvectors of 
the Hamiltonian by finding and diagonalizing optimized subspaces S. This 
process however has an intrinsic error due to the omission of basis states not 
in S. The goal of stochastic error correction is to sample the contribution of 
these remaining basis vectors. In the present investigation we use the so-called 
series expansion method up to first order. This has been explained in detail in 
[|[ [I| and therefore we restrict ourselves to the review of basic formulae. 

Let Xi be an eigenvalue and be the corresponding eigenvector of H 
restricted to the subspace S (the span of the subset of basis vectors after step 

1 We have found only one two-body bound state and one three-body bound state. Thus 
in a follow-up study one could suffice with fewer states in each sector. 
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3 of the QSE algorithm). The remaining basis vectors in the full space not 
contained in S will be denoted by \Ai), and we are interested in sampling the 
contribution of these vectors. We assume that Ai and \v\) are close to the 
corresponding quantities for the full Hamiltonian H . If we now expand in 
powers of 

«^|ff|A> - Ai)" 1 (5) 
we find that the first order correction to Ai is 

dA l-~2^ (A i \H\A i )-X 1 W 
i 

and the first order correction to \v\) is 

U„ \ _ \ " \AjUAj\Hjvi) \ \^ \v ] )(v ] \H\A,){A l \H\v 1 ) (7] 
i j^l i 

The sum over all vectors \Ai) is too unwieldy to perform exactly. Instead we 
sample the sum stochastically in the following manner. Let P(^4triai(i)) be the 
probability of selecting |^4 t riai(i)) on the i th trial. Then 

ii m 1 ^i|g|A tr iai W )(A r iai(,)|gki) 
1_ ™coAf 2^ «A trial(i) |/r|A trlal(i) >-A 1 )P(A trial(i) ) y°> 

i=l,...,N 

and 

(«, \S Vl ) = - lim 1 V H^)(A^\H\ V1 ) 

i 

, ,• 1 \ - \ - (^i^-)(^|g|A t , ia i ( ,))(^t, ia i W |gki) 

Jv" oo N 2^2^ (A J -Ai)«A tPial(i) |ff|A trial(i) >-Ai)P(A trial(i) ) 

for any fixed vector |iy). 

5 Two-body bound state 

Before presenting the results of our calculations, we first summarize what we 
expect in the weak-coupling non-relativistic limit. For the two-body bound 
state we consider all two-particle Fock states with zero momentum, 

\n,-n)=aiai n \0). (10) 
These states satisfy the orthogonality and normalization conditions 

(n 1 , -n'\ n, -n) = 5 n , n i + £„,_„'. (11) 
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At lowest order in A the : <j) e : interaction is not relevant to the two-body 
bound state problem. The effective interaction Hamiltonian is therefore 



Hi = -X : (ft 4 : (12) 

= _3A \p aL„ i a t _„ 2 a„ 3 a„ 4 

4i L^l 7 „2 X 2 \ 1/3 7 „2,2\ 57*7 „2,2 \ 1/4 / „2,-2 \ 1/4 ' 

„ 1+ „ 2+ „ 3+ „ 4=0 (V+^J (>+^) (m 2 +^-J 

In the non-relativistic limit we can approximate Hj as 

H I=—$L E a-nX-n^aru- (13) 

ni+ii2+ii3+ri4=0 

It is convenient to define the position-space basis vectors 

l*> = 7fe E e <n £*|n,-n). (14) 

n=0,±l,... 

Since we are considering a periodic box, 

| or) = |x + 2L) (15) 

and the position-space vectors satisfy 

(x'\x)= [5{x' -x + 2jL) + 6{x' +x + 2jL)}. (16) 

j=o,±i,... 

For any function ^(x), we can define 

|*) = J dx^(x) |x) . (17) 
If |*) is to have unit norm, then we require 

rL 

dx\y(x)\ 2 = \. (18) 



-L 

In our position-space representation of the two-particle subspace, we note that 

J2 *W \x) = 2 |x) (19) 

n=0,±l,... 

and 

E 4^4«nl*> = -2|@). (20) 

n=0,±l,... 

Upon restricting our Hilbcrt space to the zero-momentum two-particle sub- 
space, we can simplify the interaction Hamiltonian by the substitution 

a -n 1 a -n 2 a ™3 a «4 ~ > Xy a li 1 a ^-n 1 a n 2 a -n 2 - (21) 
ni+«2+i3+n4=0 ni,ri2 



G 



We observe that 



a ni a - ni a n 2 a-n 2 \%) = 4L6(x) \x) 



The full non-relativistic Hamiltonian, 



n=0,±l 

can now be represented as a Schrodinger operator 
We now consider the eigenvalue equation 



Using the ansatz 



we find 



When 



* 2 (x) cx e~ b|2;| +ce +6 l x l, 
6 > 0, 



_ 3AL 1+c 

c = e m i-c. 

j _ 3A 1+c 
" 2jU 1-c' 

£? = -£. 



3AL 

e « 1, 



we have 



3AL 3AL 

c = e m (l + 0(e )), 

_3AL 

6=§(l + 0(e M )), 



3AL 

f^(l + 0(e f )). 



This expression for the binding energy is consistent with the derivation 
In the infinite L limit we obtain 

3A|x| 
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6 Three-body bound state 



The analysis for the three-body bound state is slightly more involved. In the 
non-relativistic weak-coupling limit we restrict our attention to three-particle 
states with zero-momentum, 



|m,n 2 ,-ni - n 2 ) = a\^a\^_ nx _ n ^ |0) . 
We again construct position-space basis states 

at .- "l 71 " 31 ! - "2^2 ■ (-n, 1 -r. 2 )ir(- 3 ; 1 -x 2 ) 

iii,ri2=0,±l,... 



(36) 

|ni,n 2 ,-ni - n 2 ) , 
(37) 



which satisfy the orthogonality and normalization conditions 

Six', - x x + (f + 2fc)L) ( 5(4 - x 2 + f L) 
+ {xi <-> x 2 <-» -xi — x 2 } 



(x' l ,x' 2 \x ll x 2 ) = 

i,fe=o,±i,... 



The periodicity of our system implies 

\xi, x 2 ) = \x\ + §£, x 2 + |L) 

= \x\ + 2L, x 2 ) = \x\ , x 2 + 2L) . 

For any function ty(xi,x 2 ), we define 



|*) = / dxidx 2 d^(xi, x 2 ) \x\, x 2 ) 
Jr 



where R is the region 

-L <xi<L 
l 

If |^) is to have unit norm then we must have 



R 



-\L <x 2 < \L. 



I dxidx 2 \%>{xi, x 2 )\ 2 - 
Jr 



We note that 



53 ^4«n I*} = - 
n=0,±l,... 

With a little work we also find 

51 aL ni a ! _ n2 a n3 a n4 \x 1 ,x 2 ) 

™l+«2+«3+»4 = 



6- 



d 2 f , d 2 f _ d 2 <S/ 
dx'i dx'i dxidx2 



(38) 
(39) 

(40) 

(41) 

(42) 
(43) 

(44) 



= 4L 53 [S(xi + 2x 2 + 2jL) + (5(2xi + x 2 + 2jL) + 5( Xl - x 2 + 2jL)] \x t ,x 2 ) . 
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For the moment we will neglect the : cf> 6 : term. Later we will show that this 
term produces only a higher order correction in A. The effective non-relativistic 
Hamiltonian, 



H 



= E 

n=0,±l,... 



3A 



E aL ni al n2 a n3 a n4 , (45) 



"■l+'n.2+"-3+Tl4=0 

can now be represented as a Schrodinger operator: 



3A 



E 



S(xi + 2x 2 + 2jL) + 5{2xi +x 2 + 2jL) 
+S(xi -x 2 + 2jL) 



(46) 



In the infinite L limit we can write the energy eigenvalue equation as 



. JL ( ° 2 y + d 2 'f _ d 2 ^ 
I dx\ dx\ dxidx; 



3A 
M 2 



6(x l +2x 2 ) + S(2x 1 + x 2 ) 
+S(xi - x 2 ) 



Using the ansatz 



find 



Xj/ ^ e -b(\x 1 +2x 2 \ + \2x 1 +X2\ + \x 1 -x 2 \) 



6= 3A 



E 



9A 2 

7?" 



^From the normalization condition, 



3 A 

V3A „-^(\xi+2x 2 \ + \'2xi+X2\ + \xi-X2\) 



(47) 

(48) 

(49) 
(50) 

(51) 



We now return to the : </> 6 : term. This term produces an interaction 
proportional to 



t6(x 1 )S(x 2 ). 



^From (51) we see that 



(V a \5(x 1 )5(x 2 )\y 3 } = 0(\ 2 ), 



(52) 



(53) 



and therefore the : <fi 6 : interaction produces an effect of order 0(A 3 ). We expect 
however that the repulsive delta function will suppress ^3 near x\ = x 2 = 0. 
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7 Results 



In Table 1 we show for several values of A the mass of the one-particle state, 
771 1, and the two-body bound state, m^, which is the relativistic state that 
is continously connected to the non-relativistic two-body bound state. For 
comparison we show the data for m\ and mi from the lattice diagonalization 
study of M (BD), non-relativistic Schrodinger results derived above for infinite 
L (S) and L — 2.57r/j~ 1 (L), and Gaussian effective potential results from |^|, || 
(GEP). The remaining column of data comes from an exact diagonalization of 
the Hamiltonian restricted to two-particle Fock states (F). In order to obtain 
the ratio 7712/7711 in this approximation we take mi = /jj^] 

Tabic 1 



a/m 2 


mi /n 


BD 


7712/7711 


S 


L 


F 


BD 


GEP 


0.1 


0.995(1) 


0.994 


1.970(2) 


1.978 


1.967 


1.975 


1.916 


1.982 


0.2 


0.975(5) 


0.972 


1.92(1) 


1.910 


1.907 


1.934 


1.821 


1.938 


0.3 


0.94(1) 


0.928 


1.85(1) 


1.798 


1.797 


1.877 


1.718 


1.878 


0.4 


0.88(2) 


0.850 


1.76(4) 


1.640 


1.640 


1.808 


1.613 


1.807 


0.5 


0.77(4) 


0.720 


1.66(8) 


1.438 


1.438 


1.730 


1.520 


1.728 


0.6 


0.61(4) 




1.71(10) 


1.190 


1.190 


1.646 




1.642 



The error bars in Table 1 and throughout our discussion include estimated 
errors from Monte Carlo errors, higher order contributions in the stochastic error 
correction series expansion, and residual dependence on the ultraviolet cutoff 
A. Of these errors we find the higher order contributions in the expansion to 
be the largest source of error. Any follow-up study of this system should push 
the calculation to higher order or make use of the stochastic Lanczos method 
introduced in ||. 

The results from (BD) were carried out on a much smaller system. While 
their results for mi appear in agreement with ours, they have noted that their 
calculation for 777.2 is considerably affected by the small system size. At small 
coupling at least our own finite size effects appear to be small, as is suggested by 
the agreement of the Schrodinger results S and L. For this range of couplings the 
GEP predictions are equivalent to the two-particle Fock state approximation (F) 
in the limit of infinite L and A. This is because at these couplings 0, = mi = /i, 
where Cl is the optimal GEP variational mass parameter as defined in ||. 

As expected, our results agree with the Schrodinger results at small coupling. 
They also appear reasonably close to the two-particle Fock approximation, F, 
for all values of coupling. In our next calculation we measure how much of 
the bound state does in fact reside in the two-particle sector. Let I7712) be the 
normalized bound state and P2 be the projection operator onto the two-particle 
subspace. We define N2 to be the norm of the two-particle projection 



N 2 = \\P2\m2) 



(54) 



2 This is at least consistent since in the two-particle Fock approximation with normal- 
ordered interactions, the constituents of the bound state have no self-energy contributions. 
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In Table 2 we show N 2 for the different couplings. 



Tabic 2 



a/m 2 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


N 2 


0.9977(1) 


0.989(1) 


0.971(2) 


0.936(8) 


0.87(2) 


0.79(2) 



We can pursue the question of the bound state wavefunction further. Let 



*a(*) = ife (%\ m 2 ) , (55) 

where \x) is the two- particle position state defined in (|l4|). In Figures 1-3 we 
compare the wavefunctions ^2 with the corresponding approximate wavefunc- 
tions *|, *f for coupling values A = 0.1/i 2 , 0.2/j 2 , 0.4/r 2 . 

In our analysis we did not consider bound states beyond A = 0.6/i 2 . This 
is because our diagonalization/Monte Carlo results indicate a first order phase 
transition near A = 0.73/i 2 , where <f> develops a vacuum expectation value 

(0) = 1.02(2). (56) 

As a matter of fact, there appears to be no bound state in the broken symmetry 
phase, and the two-body and three-body bound states in the symmetric phase 
are stable all the way up to the transition point. The phase transition can be 
seen quite clearly in the behavior of the energy levels. In Table 3 we show 
the lowest two energy levels Eq and E% in the even <j> sector. We see what 
appears to be a metastable vacuum at energy E 2 for A < 0.73/i 2 becoming the 
true vacuum at energy Eq for A > 0.73/i 2 .f^] 



Tabic 3 



A//i 2 


E a 


E 2 


0.70 


-0.37(5) 


0.09(10) 


0.71 


-0.41(5) 


-0.15(15) 


0.72 


-0.40(5) 


-0.29(20) 


0.73 


-0.49(20) 


-0.42(5) 


0.74 


-0.70(20) 


-0.42(5) 



In Table 4 we show for several values of A the mass of the three-body bound 
state, 7713, which is the relativistic state that is continously connected to the 
non-relativistic three-body bound state. We were not able to find any pre- 
vious studies of this state in the literature. For comparison we show data 
from the Schrodingcr equation at infinite L (S) and two different Fock space 
approximations, Fi and F2. Fi corresponds with an exact diagonalization of 
the Hamiltonian restricted to three-particle Fock states, and F2 corresponds to 
the same thing, except with the : </> 6 : interaction turned off.^J To obtain the 

3 Since we are at finite volume the levels Eq and E2 come close but never actually become 
degenerate. 

4 This is an artificial approximation. Without the restriction to three-particle Fock states, 
turning off the : 6 : interaction would make the theory unbounded below. 
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ratio ms/mi in the two Fock approximations we again take mi = /i. 



Tabic 4 





ms/mi 


S 


Fi 


F 2 


0.1 


2.931(2) 


2.910 


2.940 


2.919 


0.2 


2.83(1) 


2.640 


2.862 


2.770 


0.3 


2.72(2) 


2.190 


2.769 


2.568 


0.4 


2.59(6) 


1.560 


2.667 


2.336 


0.5 


2.54(8) 


0.75 


2.558 


2.084 


0.6 


2.62(15) 


< 


2.444 


1.818 



We notice that the Schrodingcr approximation (S) is not quite as accurate for 
given A as it was for the two-body bound state. Nevertheless the Fock approxi- 
mation Fi still appears to follow the actual data relatively well. We observe that 
7713 is not only below the 3toi continuum threshold but also below the mi + mi 
threshold. The corresponding state therefore has a sensible interpretation as a 
stable particle. 

In analogy with before, we let \m^) be the normalized three-body bound 
state and P3 be the projection operator onto the three-particle subspace. We 
define iV 3 to be the norm of the three-particle projection. 

JV 3 = ||f%|ro 3 >||. (57) 
In Table 5 we show N3 for different couplings. 



Tabic 5 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


N 3 


0.9966(2) 


0.981(2) 


0.950(5) 


0.89(1) 


0.80(2) 


0.66(4) 



Let 



^3{xi,x 2 ) = ± (xi,x 2 \ m 3 ) . (58) 

In Figures 4-11 we compare the wavefunction ^3 with the corresponding ap- 
proximate wavefunctions \Pf, "J^ 1 , ^f 2 for A = 0.1/x 2 and 0.4^, 2 . For visual 
clarity we plot the wavefunctions over three copies of the fundamental region 
R. We use the variables 

Vi = 3EZi, (59) 
Vi = # (*! + 2x 2) , (60) 

to make clear the hexagon symmetry. To guide the eye we have drawn dotted 
lines along the points 

x\ = x 2 , (61) 
x\ — —x\ — x 2 , (62) 
x 2 = —x\ — x 2 . (63) 
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These corresponds with configurations where two of the three constituents lie 
at the same point. 

We notice that for both ^3 and ^f 1 , the inclusion of the repulsive : 6 : 
interaction produces a marked suppression of the wavefunction near the origin. 
By independently changing the coefficients of : <fi 6 : and : 4 :, we have a simple 
laboratory to study the competition between two-body and three-body forces 
in three-body bound states, a topic of some relevance and recent interest in the 
triton system. jl0| 



8 Summary 

We have considered the relativistic bound states of 4> e — cf> 4 theory in 1 + 1 dimen- 
sions. Using a recently proposed diagonalization/Monte Carlo computational 
scheme, we calculated the two-body and three-body bound state energies and 
wavefunctions. The initial diagonalization was performed using the quasi-sparse 
eigenvector method, and the subsequent stochastic error correction was calcu- 
lated to first-order using the series method. We made detailed comparisons of 
our numerical results with the available literature and several approximations, 
including non-relativistic weak-coupling solutions and various truncated Fock 
space approximations. 

Applications to more complicated higher dimensional theories can be done 
in the same manner but require larger QSE subspaces and going beyond the 
first-order term in the series method or use of the stochastic Lanczos method. 
But even from the simple examples presented here, the potential for studying 
relativistic bound states with diagonalization/Monte Carlo methods is apparent. 
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Three-body wavefunction * g at Wfi 2 = 0.1 




Figure 4: Three-body wavefunction ^ 3 at A = O.lfi 2 . 
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Figure 5: Approximate three-body wavefunction \&| at A = O.l^ 2 . 
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Three-body wavefunction at A/|i 2 = 0.1 
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Three-body wavefunction * 3 at A/^ 2 = 0.4 




Figure 8: Three-body wavefunction ^3 at A = 0.4/x 2 . 
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Figure 9: Approximate three-body wavefunction <i?f at A = 0.4/x 2 . 
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Three-body wavefunction *p|li at A/|i 2 = 0.4 
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